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It has been known for a long time that light carries both linear and angular momenta parallel to 
the direction of propagation. However, only recently it has been pointed out that beams of light, 
under certain conditions, may exhibit a transverse spin angular momentum perpendicular to the 
propagation direction. When this happens, the electric field transported by the light rotates around 
an axis transverse to the beam path. Such kind of fields, although deceptively elusive, are almost 
ubiquitous in optics as they manifests in strongly focused beams, plasmonic fields and evanescent 
waves. In this work we present a general formalism describing all these phenomena. In particular, 
we demonstrate how to mathematically generate a wave field possessing transverse spin angular 
momentum density, from any arbitrarily given scalar wave field, either propagating and evanescent. 


INTRODUCTION 

The electric field of a monochromatic electromagnetic 
wave of angular frequency uj spans, at any fixed point in 
space, an elliptical trajectory in a period T = 2tt/uj [I]. 
In the points where the ellipse reduces to a circle, the 
light is said to be circularly polarized. When the light 
propagates in the direction perpendicular to the plane of 
such circles, the uniform motion of the electric field vector 
generates a longitudinal spin angular momentum, loosely 
dubbed “helicity” ||2j. However, there are instances where 
the light is circularly polarized in a plane containing the 
direction of propagation. In these cases the light carries a 
transverse spin angular momentum density mg. When 
the light wave exhibits both transverse and longitudinal 
electric field components, the transverse spin angular mo¬ 
mentum can be either dependent or independent of the 
longitudinal one m- 

In this work we aim at establishing a perfectly gen¬ 
eral theory of wave fields displaying transverse spin angu¬ 
lar momentum (AM). We begin with studying particular 
solutions of Maxwell’s equations possessing both trans¬ 
verse and longitudinal components of the electric field 
vector. The kinematics of the latter is investigated by 
means of suitably defined Stokes parameters pHL This 
parametrization permits us to establish general sufficient 
conditions for the existence of transverse spin AM. Un¬ 
der certain circumstances, such conditions coincide with 
the Cauchy-Riemann equations satisfied by holomorphic 
complex functions m • This novel result establishes 
an intriguing and previously unnoticed connection be¬ 
tween the theory of complex functions and the transverse 
spin AM of light. Finally, to illustrate our findings, we 
work out two detailed examples involving non-diffracting 
Bessel beams and evanescent waves. Caveat : In the re¬ 
mainder, we adopt units where the electric field has the 
dimensions of a wave-number, namely the inverse of a 
length. 
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TRANSVERSE SPIN ANGULAR MOMENTUM 

In the Coulomb gauge m the electromagnetic vector 
potential A(r, t) is purely transverse, namely V • A = 
0, and satisfies the wave equation DA = 0. A plane- 
wave solution of these equations takes the form A(r, t) = 
n exp(ik • r — iuot), where k is the either real- or complex¬ 
valued wave vector ma. uj = c|k| = ck is the angular 
frequency and n is a constant three-vectors such that 
k • n = 0. In the Cartesian reference frame (x, y , z) the 
latter condition can be written as 

k • n — k x n x H~ kyny -1- k z n z — 0. (1) 

Apart from the trivial solution n = (0,0,0), Eq. 0 
also admits three elementary solutions of the form ni = 
(—ky, k x , 0), n 2 = (-k z ,0,k x ),n 3 = (0 ,-k z ,k y ), where 
one of the three Cartesian components of n is chosen to 
be zero. The first solution ni generates a plane-wave 
mode whose electric field is purely perpendicular to the 
z -axis H23- These modes are commonplace in the theory 
of electromagnetic waveguides m • For instance, since 
the condition k • ni = 0 is fulfilled irrespective of the 
value of k z , one can choose k = ki = (k x ,k y , 0 ) and 
obtain the so-called transverse electric (TE) modes of a 
rectangular waveguide 0. The remaining two solutions 
112,113 generates plane-wave modes with both longitudi¬ 
nal ( z ) and transverse (either x or y, respectively) compo¬ 
nents of the electric field. As an example, when choosing 
k = k 2 = (fc x , 0, k z ) and k = k 3 = (0, k y , k z ) with n = n 2 
and n = n 3 , respectively, one attains the transverse mag¬ 
netic (TM) modes of a rectangular waveguide. 

Consider now n = n 3 . In this case the vector potential 
can be written as 

A(r, t) = (0, —k z , k y ) exp(ik • r — iuot) 

= — i(0, — d z , d y ) exp(zk • r — iut). (2) 

Let ^(r, t) be an arbitrary solution of the wave equation, 
namely □V’ = 0. By definition, this field can always be 
expressed as a superposition of plane waves. Therefore, 
using Eq. 0 one can write A = (0, — dz'i/j , d y i/j), with A 
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automatically satisfying both V • A = 0 and DA = 0. 
The electric and magnetic fields generated by A are 

E = — A = (0, d z ip, —dytp ), (3) 

B = V x A = (dyip + d 2 z xf>, -d x d y ^, -dM ), (4) 

where, A = 9 t A and ip = dt'ip- It is interesting 
to note that the electromagnetic fields above can also 
be obtained from the Hertz magnetic vector potential 
n m = (-i/j/no, 0,0), with A = yitoV x II m dHUIS]. 

When pj represents a time-harmonic wave of the form 
^(r, t) = U(r) exp (— iuot) /u, the physical electric field 
vector 

E jR (r, t) = Re [—i(0, d z U, — d y U) exp (— iut)] 

= p(r) cos (cat) + q(r) sin (jut) (5) 

sweeps out an ellipse lying in the yz- plane, where p = 
(0, d z v , —d y v) and q = (0, —d z u, d y u ), with u(r) and v(r) 
denoting, respectively, the real and the imaginary part of 
U = u + iv [I|. When the latter is given, the magnitude 
and the orientation of the axes of the ellipse at a fixed 
point r = (x, ?/, z) in space, are determined by the coor¬ 
dinates x, y and 2 of the point. In particular, if at r = 1*0 
p(i*o) • q(ro) = 0 and p 2 (i*o) — q 2 (ro) = 0 are fulfilled, 
then the wave is circularly polarized in the yz -plane [l] 
and the spin AM density is purely transverse at r = r 0 . 
In other points on the yz- plane, the polarization can be 
completely characterized by the spatially-varying Stokes 
parameters So(r), Si(r), 62 ( 1 *), Ss(r) defined in terms of 
the vectors p, q, as 

50 = p 2 z + qI + pI + qI 

= (< d y v) 2 + ( d y u) 2 + ( 1 d z v) 2 + ( d z u ) 2 , (6) 

51 = p\ + q 2 z - P 2 V ~ ql 

= ( d y v ) 2 + ( d y u ) 2 - ( d z v ) 2 - ( d z u ) 2 , (7) 

5 2 = 2 (PzPy + q z q y ) 

= — 2 ( d z ud y u + d z vd y v ), (8) 

5 3 = 2(p z q y — p y q z ) 

= 2 ( d z ud y v — d y ud z v ), (9) 


where we have used the notation by Dennis [[16 . 

The first constraint p(ro)-q(i*o) = 0 for the existence of 
circular polarization at r = r 0 , can be expressed in terms 
of u , v as: d z ud z v + d y ud y v = 0 . A possible solution of 
this equation is 


(d z v) (r 0 ) = (~d y u) (r 0 ), (d y v) (r 0 ) = ( d z u ) (r 0 ), (10) 


where the notation emphasizes the fact that the equa¬ 
tions in ( 10 ) are supposed to become equalities only in 


the isolated point ro, and not necessarily everywhere on 
the yz- plane. Then at r = r 0 , p = —(0, d y u , d z u ) and the 
second condition p 2 — q 2 = 0 is automatically satisfied. 
Moreover, from Eqs. ([^[9| it follows that Si = S 2 = 0 
and S 3 /So = 1 , as expected for a circularly polarized 


wave. 


Equation (10) is the first main result of this work and 
establishes a necessary and sufficient condition for the 
existence of transverse spin AM density in fields of the 
form ffl- This concept can be further extended by tak¬ 
ing for V = u+iv, a holomorphic function of the variables 
y and z. Then, the Cauchy-Riemann equations satisfied 
by V to be holomorphic, read d z v = — d y u and d y v = d z u 
HH> which coincides with Eq. ( |10| ) at r = ro- Therefore, 
when U is holomorphic the light carries a transversely 
spinning electric field (or, synonymously, transverse cir¬ 
cular polarization) uniformly distributed everywhere on 
the yz- plane. Elementary solutions of the monochro¬ 
matic Helmholtz equation (V 2 + k 2 )U = 0, which are 
also holomorphic with respect to the variables y and z, 
must have the form U(r) = H(y, z) exp (ikx), where H is 
an arbitrary harmonic function on the yz- plane, namely 
(<9 2 + dl)H = 0 . Of course, by definition of harmonic 
functions, H cannot be bounded on the yz -plane and, 
therefore, it does not admit a finite norm. However, this 
does not prevent the existence of finite functions approx¬ 
imating H (see, e.g., discussions in mini). 


A simple, practical “recipe” to mathematically obtain 
uniform circular polarization everywhere on the yz- plane, 
can be given by simply noting that if /(£) is an arbitrary 
smooth real function with respect to the real variable 
£, then f{z + iy) = a(z, y) + ib{z , y) is by definition a 
holomorphic function with respect to the variables y and 
z m- Now, with the harmonic real functions a and 6 , 
one can build the two reciprocally orthogonal vectors p = 
(0, a, b) cos(kx) and q = (0, — 6 , a) cos(kx). Substituting 
these vectors in Eq. ([5|, one obtains an electric field with 
uniform transverse circular polarization on the yz- plane. 

The physical content of our theory can be further elu¬ 
cidated by calculating the quantity • E^, key to the 
magnetic circular dichroism manifested by molecules in 
DC magnetic fields HU From Eq. © it follows that 
E R • E^ = uj (p x q). Comparing this result with Eq. © 
we immediately find 


E«.E fl = -|(5 3 ,0,0). (11) 

This means that our wave field © couples to molecules 
subjected to a DC magnetic field in the x direction and 
that the coupling is maximum (S 3 = ± 1 ) at the points 
where the electric field vector sweeps out a circle on the 
yz- plane. 

The main virtue of the fields of the form <H§ re¬ 
sides in their straightforward physical significance that 
leads to the establishment of an intriguing connection be¬ 
tween the distribution of transverse spin AM density in 
light waves and the theory of complex functions. How¬ 
ever, more general wave fields with nonzero longitudi¬ 
nal electric field vector component can be obtained by 
considering an electromagnetic vector potential of the 
form A = A 2 + A 3 , where A 2 = (—d z (p, 0,3^0) and 
A 3 = ( 0 , —d z ip,d y ip), with 0 (r,£),^(r, t) being two in¬ 
dependent solutions of the wave equation. This choice 
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yields the electric field vector 


E = (d z <p, d z ip, -d x 4> - dyip), 


( 12 ) 


which simultaneously displays both nonzero transverse 
and longitudinal Cartesian components. The expres¬ 
sion in Eq. (12) is interesting because it also implies 


that a wave of the form ( 12 ) has necessarily a vanish¬ 


ing longitudinal component when the light is uniformly 
circularly polarized everywhere in the xy- plane. This 
can be seen again by using complex calculus. Indeed, 
if we choose = icj) with 0(r, t) = U(r) exp (— iut) /cj, 
then we have p = (d z v, d z u , — d y u — d x v ) and q = 
(— d z u, d z v, d x u — d y v). The conditions p • q = 0 = 
p 2 —q 2 to be fulfilled in order to have circular polarization 
in the plane spanned by p and q, now read d y v = d x u, 
and d x v = —d y u. The latter relations coincide with the 
Cauchy-Riemann equations for holomorphic functions in 
the variables x and y. When satisfied, these equations 
lead to p z = 0 = q z , namely to a vanishing longitudinal 
electric field vector component. This case illustrates once 
more the strict connection existing between transverse 
circular polarization of light and the theory of complex 
functions. 


The temporal evolution of E jR (r, 0, t) as a function of the 
scaled coordinates z/zo and y/ro, superimposed on a den¬ 
sity plot of the norm of the field (E^ • E^) 1 / 2 , is shown 
in Fig. [l] below. In the remainder, x and y coordinates 
are normalized to ro = a/k r , with a « 2.405 being the 
smaller solution of the equation Jo (a) = 0 , while the z 
coordinate is normalized to zo = l/k z . Moreover, in all 
plots we choose i9q = 15°. 


6Jt=0 a> t = ZZ/2 (Ot = ZT 



- 2-10 1 2 - 2-10 1 2 - 2-10 1 2 


z/z 0 z/z 0 z/z 0 


FIG. 1. Instantaneous distributions of the real-valued elec¬ 
tric field E R (r,0,t) ( [15] ). At y/ro ~ 1 the electric field spins 
in the yz- plane. Here ro = a/k r gives the central core spot 
size of the beam, with a ~ 2.405 such that Jo (a) = 0 and 
$o = 15°. For comparison, see [9]. 


APPLICATIONS OF THE FORMALISM 

The scalar field ^(r, t) considered up to now, was com¬ 
pletely arbitrary. This choice permitted us to find some 
perfectly general fundamental relations. However, for the 
sake of definiteness, we now specialize in the study of a 
scalar zeroth-order Bessel field and of an evanescent plane 
wave. 


A. Bessel field 


Consider the scalar Bessel field 

^(r, t) = J 0 (rk r ) exp (izk z ) exp ( ? ( 13 ) 


where k r = k sin $ 0 and k z = k cos$o, with $ 0 denoting 
the aperture of the Bessel cone and x = rcos0, y = 
r sin 0 defining the polar coordinates (r, 6 ) in the xy- 
plane m- From Eq. © it follows that U = 
Jo {rkr) exp (izk z ) with u = Jo (rk r ) cos (zk z ) and v = 
Jo (rk r ) sin (zk z ). The conditions for circular polariza¬ 
tion Eq. ([IQ]) reduce therefore to the single equation 


cot'd 


y J\ (rk r ) 
r J 0 {rk r ) ’ 


(14) 


which can be numerically inverted to find the points ro = 
(xo,yo,z) where Eq. (14) becomes an identity. 

The real electric field generated by using Eq. 

|3| , can be expressed in polar coordinates (r, 6 ) as 


E^(r, 0,t) = y k z Jo ( rk r ) cos (zk z — ut) 

+ z k r sin 6 Ji (rk r ) s'm (zk z — wt). (15) 


The polarization distribution on the yz -plane is char¬ 
acterized by the Stokes parameters 


So — <y 2 + /3 2 , 

(16) 

S\ = a 2 — (3 2 , 

(17) 

s 2 = o, 

(18) 

S 3 = 2a/3, 

(19) 


where a = (y/r) k r Ji(rk r ) and /3 = k z Jo(rk r ). Equation 
(18) shows that diagonal polarization never occurs in this 
kind of beams. The occurrence of transverse circular po¬ 
larization is naturally marked by the conditions S 1 = 0 
and S 3 = So- These are clearly satisfied when a = /?, 
namely when Eq. (fl4| becomes an identity. The ratio 
S 3 /So is plotted as a function of the scaled coordinates 
x/ro and y/ro in Fig. [ 2 ] This figure shows that one has 
S 3 /So = ±1 in the tail of the central core spot of the 
beam, namely for x 2 + y 2 ~ Tq. Moreover, a discontinu¬ 
ity occurs around x/ro = ±1 and y —> 0 , as illustrated in 
detail in Fig. [3] 

Despite of the unusual spatial distribution of S 3 /So, 
the projection of the beam upon the xy -plane displays the 
typical Bessel form. This can be seen by calculating the 
time-averaged energy density w(x,y), which is conserved 
along 2 ) and is given by 


™0r,y) = !(E*-E + c 2 B*-B), (20) 

where eo is the vacuum permittivity [10]. The plot of 
w(x,y) normalized to the central value wo = w( 0 , 0 ) is 
shown in Fig. [4] 

Another relevant quantity that permits to characterize 
the beam, is the distribution of the linear momentum P 
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FIG. 2. Ratio S3/So as a function of the scaled coordinates 
x/ro and y/ro. The electric field vector manifests a (local) 
transverse circular polarization at the points where S3/So = 
± 1 . 


FIG. 4. Time-averaged energy density w(x,y)/wo as given 


by Eq. (20), with 7*0 = a/k r and $0 = 15°. 


of the field, which naturally splits into a canonical and a 
spin part asP = Pc + Ps, where 

P c = Irn [E* • (V) E + c 2 B* • (V) Bl , (21) 

4CJ 


Ps = |v x S, (22) 

with 

S = — Im (E* x E + c 2 B* x B) 

4 cj v ’ 

= S e + S#, (23) 

being the spin AM density [20] . A straightforward cal¬ 
culation shows that P c = z P Ez , P s = z P Bz and 
S = -kS x + y S y . This implies that S • P = 0, namely 
the helicity of the beam is identically zero: 

h = —— Im(E* -B) =0. (24) 


This peculiar phenomenon also occurs for surface- 
polariton fields and focused beams mm- The spatial 

1.0 


0.5 


o 


- 0.5 


- 1.0 



FIG. 3. Expanded 3D view of the discontinuity occurring in 
S 3 /So around x/ro = 1 and y ~ 0. 


patterns of P^ and P s are shown in Fig. [5] It is of worth 
noticing that Ps z takes negative values in some annular 
regions of the xy- plane, but the sum Pq z + Psz > 0 is 
nonnegative everywhere. This peculiar phenomenon is 
well known for the transverse spin AM density in evanes¬ 
cent waves (see [6: and discussion later). 



FIG. 5. Distributions of the ^-component of the time- 
averaged canonical (left) and spin (right) line ar mo menta, 
Pc and P s , respectively, as given by Eqs. ( 2l][22 ). The 
occurrence of negative values for Ps z , should be noticed. 


Although the helicity (24) of the wave is identically 
zero, the spin AM density is not. Since S = x + y S y , 
it is clear that the spin AM for the wave (15) is purely 
transverse. A straightforward calculation shows that 
S E = x SEx and S B = xS Bx + y S By , with \S Bx \ < 
\S By \ ~ \S Ex \- Figure [6] shows the vector field distribu¬ 
tion on the xy- plane of (from left to right) S E , S B and 
S = Se + S B , respectively, superimposed on a density 
plot of the norm of the corresponding fields. Both the 
electric and magnetic fields contribute to the transverse 
spin AM density, the two terms being in quadrature. As 
we shall show soon, the situation is drastically different 
from the one occurring with evanescent waves, where only 
the electric field contribution is nonzero. 
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FIG. 6. (Left to right) Distributions of the the spin AM den¬ 
sities vector fields S e, S b and S = Se + Sb on the xy- plane, 
calculated from Eq. ( [23] ) . All plotted fields are normalized to 

Wq/u. 


FIG. 7. Instantaneous distributions of the real-valued elec¬ 


tric field E K (y,z,t) as given in Eq. (28) and evaluated for 


£ = 7r/12. In each point the electric fie" 
wise. 


d vector rotates clock- 


B. Evanescent waves 


Consider now the evanescent field of an inhomogeneous 
plane wave exponentially decaying in the positive z di¬ 
rection and propagating along the y axis m : 


</>(r, t) = exp (—zk sinh ( + iyk cosh () 


n —iu;t 


W 


(25) 


where the real parameter ( fixes the scale of inhomogene¬ 
ity. For C — 0 the wave is purely propagating. The elec¬ 
tric and magnetic fields obtained by substituting 0(r, t) 
into Eqs. i@, are given by 


E = — ck 2 <p (—iy sinh ( + z cosh (), 
B = - k 2 c 


>x. 


(26) 

(27) 


The corresponding real electric field is written as 
E R (y,z,t) = — k exp (—zk sinh 

y sin (yk cosh ( — wt) sinh ( 

+z cos (yk cosh ( — ujt) cosh ( . (28) 

Figure [ 7 ] shows the temporal evolution of E R (y,z,t) as 
a function of the scaled coordinates ky and kz, superim¬ 
posed on a density plot of the norm of the field. 

For an evanescent wave, the energy density depends 
on 2 solely and for z > 0 it decays as w(z)/wo = 
exp (—2zk sinh (). After a straightforward calculation, 
from Eq. (23) we obtain Sb = 0 and Se = S, where 


UJ 


— S = x exp (— 2 zk sinh Q tanh (. 
w 0 


(29) 


This expression (apart from a different choice of the 
transverse axis) is in agreement with [6]. The canon¬ 
ical and the spin momenta can be calculated by using 
Eqs. ( |2lp2| and the result is 


—P c = y exp (—2 zk sinh C) cosh C, 
wo 


(30) 


—P s= — y exp (—2z k sinh 0 sinh ("tanh (. (31) 

w 0 


Also these results are in agreement with Eqs. (8-9) in [§]. 
It should be noted that the “backward” momentum P s 
is compensated by the forward momentum P c and their 
sum is eventually forward: 


— (Pc + P s) = y exp (—2 zk sinh Q sech (32) 

w 0 

Finally, the Stokes parameters are written as 

50 = k 2 cosh(2C) exp (—2 zk sinh (), (33) 

51 = k 2 exp (—2 zk sinh (), (34) 

^2 = 0, (35) 

S 3 = — k 2 sinh(2C) exp (—2 zk sinh (), (36) 


with S 3 /So = — tanh(2£). This means that for ( > 1 the 
electric field has (almost) transverse circular polarization 
uniformly over the yz- plane. 

More complex evanescent waves may be generated by 
interference. For example, suppose to illuminate the ver¬ 
tical edge of a parallelepiped immersed in a liquid with 
higher refractive index. Under suitable conditions, one 
obtains inside the solid two perpendicular evanescent 
waves, namely </>(r, t) given in Eq. (25) and <p±(i *,£) = 
exp (—yk sinh ( + izk cosh () exp(—The time- 
evolution of the interference pattern generated by sum¬ 
ming and </>_l is shown in Fig. [8j 


oj t = 0 cot = sr/4 cot = jt/2 



024680246802468 

k z k z k z 


FIG. 8. Instantaneous distributions of the real-valued 
electric held E R (y,z,t) generated by (j)(r : t) + c/)±(i?,t) with 
C = tt/12. 






























































































6 


DISCUSSION AND CONCLUSIONS 


The reason why transverse circular polarization is most 
easily observed in strongly focused beams, is evident from 
Eq. §. For a paraxial fundamental Gaussian beam with 
angular aperture < 1, one has that dzip ~ 1 + 0(0%) 
and dy^ip ~ O(0o). In this case the transverse and lon¬ 
gitudinal components of the electric field vector do not 
have a comparable magnitude. Alternatively, one may 
consider using waveguides, where TM modes with trans¬ 
verse elliptical polarization are created by superposition 
of plane waves m EH- As a matter of fact, it is suf¬ 
ficient to consider the interference of two plane waves 
only, to obtain non-paraxial optical fields with nontriv¬ 
ial polarization patterns comprising transverse circular 
polarization, as noticed in [22 and [8]. 


mentum. Of course, one may think to carry out an analy¬ 
sis similar to the one presented in this work, but studying 
instead transverse orbital AM. We are currently investi¬ 
gating along this direction. 

In conclusion, we have presented a perfectly general 
theory of light carrying an electric vector field circularly 
polarized in a plane containing the main axis of propa¬ 
gation. The novelty of our approach resides in its “uni¬ 
versal” character, which provides for a unifying view of 
seemingly different wave propagation phenomena. The 
success of such unification is made manifest in the two ex¬ 
amples reported in this paper, where it is shown that for 
both propagating (Bessel) and evanescent (plane wave) 
fields, the underlying mechanism generating transverse 
spin AM is the same. Last but not least, our treatment 
reveals a somewhat hidden connection between light with 
circular polarization and the theory of complex functions. 


For the Bessel field (13) we have found that the or¬ 
bital linear momentum is purely longitudinal, namely 
P c = z Pcz- This implies that the orbital part of the an¬ 
gular momentum density is purely radial: Jo = r xPc = 
(?/x — xy)Pcz leading to a null total orbital angular mo¬ 
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